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ABSTRACT 

Radio, infrared, and optical observations of the 2006 eruption of the symbiotic recurrent nova RS 
Ophiuchi (RS Oph) showed that the explosion produced non-spherical eject a. Some of this ejected 
material was in the form of bipolar jets to the east and west of the central source. Here we describe X- 
ray observations taken with the Chandra X-ray Observatory one and a half years after the beginning 
of the outburst that reveal narrow, extended structure with a position angle of approximately 300 
degrees (east of north). Although the orientation of the extended feature in the X-ray image is 
consistent with the readout direction of the CCD detector, extensive testing suggests that the feature 
is not an artifact. Assuming it is not an instrumental effect, the extended X-ray structure shows hot 
plasma stretching more than 1,900 AU from the central binary (taking a distance of 1.6 kpc). The 
X-ray emission is elongated in the northwest direction — in line with the extended infrared emission 
and some minor features in the published radio image. It is less consistent with the orientation of 
the radio jets and the main bipolar optical structure. Most of the photons in the extended X-ray 
structure have energies of less than 0.8 keV. If the extended X-ray feature was produced when the 
nova explosion occurred, then its 1" .2 length as of 2007 August implies that it expanded at an average 
rate of more than 2 mas d _1 , which corresponds to a flow speed of greater than 6,000 km/s (d/1.6 
kpc) in the plane of the sky. This expansion rate is similar to the earliest measured expansion rates 
for the radio jets. 

Subject headings: Binaries: symbiotic — stars: individual (RS Ophiuchi) — novae, cataclysmic vari- 
ables — X-rays: binaries 



1. INTRODUCTION 

Novae are the most common major stellar explosions 
in the universe. The basic properties of these explo- 
sions are well explained by models in which a thermonu- 
clear runaway at the base of the envelope of an ac- 
creting white dwarf (WD) causes the WD to brighten 
to near the Eddington luminosity and eject its en- 
velope in a spherical outflow. An increasing num- 
ber of observations, however, are indicating that the 
eject a from many novae are anything bu t spherical (e.g., 
ISlavin et alJUOQa IGill fc O'Brienl l2000l : iKrautter et all 
120021 : lHarman fc ? Brie3l2QQ3[ ). ev en in the X-ray s (de- 
tected a century after outburst bv iBalmanl l2QQ5h . The 
properties of the outflows triggered by nova explosions 
on massive WDs are particularly important because they 
may influence or reveal the amount of mass lost per ex- 
plosion, and thereby the ability of the stars to become 
type la super novae. 

Since the recurrence time for a nova explosion in a 
given system is inver sely related to the mass of the 
WD (jYaron et al.l 120051 ). interacting binaries that con- 
tain WDs with masses greater than approximately 1.1 
M can have recurrence times of less than 100 years; 
these systems are the so-called recurrent novae. In sym- 
biotic recurrent novae - of which only 4 are known: RS 
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Oph , T CrB, V3890 Sgr and V745 Sco - the accreting 
WD is embedded in a dense nebula fed by the wind of 
the red-giant companion. RS Oph went into outburst 
for the sixth recorded time (previous outbursts occurred 
in 1898, 1933 1958 1967, 1985) on February 12th 2006 
(N arumi et a l. 2006). The astronomical community ob- 
served this outburst at virtually every wavelength; the X- 
ray range was extremely well cov ered, with observations 
bv R XTE dSokoloski et al.l l2006h. Swift (|Hachisu et all 
120071: iBode et alJl2006l)~/mn^ra, and XMM-Newton 
(|Ness et al.l 120071 : iNelson et al.ll2008[ ). 

Observations of RS Oph indicate that the morphol- 
ogy of the expanding shell is complex and highly asym- 
metric. Within a few weeks of the start of the 2006 
eruption, X-ray spectra showed that the blast- wave evo- 
lution deviated from self-similar , sp herical expansion 
(jSokoloski et all 120061 : [Bode et al.ll2006f ). Bipolar struc- 
ture along the east-west direction was clearly seen in 
radio images taken a few weeks into the 2 006 out- 
burst (|Q'Brien et all 120061 : iRupen et all |2QQ8[ ): several 
months after the start of the outburst, observations 
showed that the radio jets were composed of highly 
collim ated outflows that p owered synchrotron-emitting 
lobes (Sokol oski et aLll2Q08l ) . At optical wavelengths, [Ne 
VJA3426 images 5 months after the start of the outburst 
revealed an exp anding nebular re mnant with a double- 
ring structure (|Bode et al.l 120071 ). The major axis of 
the double ring was also oriented east- west. Near in- 
frared (NIR) interferometric observations also revealed 
extended asymmetric structure in the early expanding 
shell — but along a differe nt direction. On day 5 .5 after 
the start of the outburst, iChesneau et al.l (|2007l ) found 
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that the NIR emission arose from a rapidly expanding 
elliptical region with major axis oriented from northwest 
to southeast. NIR obser vations in the first month of the 
outburst by lLane et all (|2QQ7l ) also revealed an asym- 
metric structure with a position angle (P. A.) of approx- 
imately 120°/300°. 

In this letter, we describe the first detection of ex- 
tended X-ray structure from RS Oph. Taken one and half 
years after the start of the 2006 eruption, our Chandra 
observations reveal a narrow X-ray feature along the di- 
rection of the NIR extended emission. We describe the 
observation in section £[2j detail the analysis and results 
in section §3j and discuss possible interpretations of the 
results in section [JU 

2. OBSERVATIONS 

On 2007 August 4, the Chandra X-ray Observatory 
(CXO) observed RS Oph for 90.1 ks using the ACIS- 
S S3 back illuminated chip with a 1/4 subarray (Obsld 
7457, start time 06:11:26 UT), in Timed Exposure mode 
and Faint telemetry format. The chip was read out every 
0.94104 s, with the time between each read consisting of 
a 0.9 s exposure and 0.04104 s for charge transfer. RS 
Oph produced an average count rate of 0.0682 counts 
s _1 , or 0.0642 counts per frame integration time. The 
configuration we used ensured a minimal pile-up fraction 
(PIMMS 5 calculates a negligible pile- up fraction of 2%). 

3. DATA ANALYSIS AND RESULTS 

3.1. X-ray spectrum 

With a source count rate of 0.0682 c s _1 (see §[2j), the 
background subtracted spectrum contained 6252 counts. 
The X-ray spectrum of RS Oph in quiescence is com- 
plex. Here we use a preliminary fit (x^=1.35 for 154 
d.o.f.) consisting of two soft, relatively unabsorbed 
thermal plasma components (mekal) and a hard, heav- 
ily absorbed cooling flow (mkcflo w) component (usin g 
the model names found in XSPEC: lArnaud et alJ ll996). 
This fit adequately models the observed spectrum for the 
purposes of the PSF deconvolution processes. We will 
present a full analysis and interpretation of the spectrum 
in a forthcoming paper (Nelson et al. in preparation). 

3.2. Extended X-ray structure 

Since observations at other wavelengths suggested that 
any extended structure would be on the order of an arc- 
sec in len gth, we applied a Subpixel Event Reposition- 
ing (SER: lLi et all 120031 . 120041 ) technique to obtain the 
best possible spatial resolution. On the basis of a phys- 
ical model of the photon-CCD interaction, SER uses 
the available information concerning photon event charge 
splitting to reconstruct the most likely entrance positions 
of photons within individual CCD pixels. This proce- 
dure potentially improves the effective Chandra/ ACIS 
pixel size from 0".5 to as small as 0".3, albeit with a 
non- uniform shape and a strong dependence o n event 
charge split (i.e., event grade); see lLi et al] (|2003l ) for de- 
tails. Due to the significant number of counts detected, 
we used a spatial bin size of 0".125 in constructing the 
X-ray images analyzed here, so as to make best use of 
the potential improvement in spatial resolution afforded 
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by SER. After applying SER, we performed a maximum 
likelihood deconvolution on the image using a synthetic 
Point Spread Function (PSF) appropriate for the source 
off-axis angle and X-ray spectrum. Altogether, the pro- 
cess for generating the final image consisted of the fol- 
lowing steps: 1) extract and fit the spectrum ; 2) u se the 
Chandra Ray Tracer fChaRT: ICarter et al.ll2003h . with 
the spectral model from #1, to trace the rays through 
the Chandra X-ray optics, and then use MARX (Model 
of AXAF Response to X-rays) to project those rays onto 
the detector and create a synthetic PSF 3) apply SER 
to the pipeline-processed event file ; 4) use the synthetic 
PSF from step #2 and the SER-processed event file from 
step #3 to perform maximum likelihood deconvolution. 

The maximum likelihood deconvolution of the 
Chandra observation reveals a feature extending north- 
west of the central binary (a=17h 50m 13.2s, £=-06° 42' 
28.2") at a P.A. (east of north) of approximately 300° 
(Fig. [1] and [2] ). The length of the extended feature, as 
measured from the centroid of the central X-ray emission 
to the location where the intensity falls below 10 counts 
per // .125 pixel, is 1.2±0.3 arcsec. We took the uncer- 
tainty on the length of the extended structure to be the 
best achievable SER pixel size, i.e., 0".3. A comparison 
among images from four different energy bands (selected 
to contain roughly 1,000 counts each) showed that the 
narrow X-ray feature is primarily composed of photons 
with energies less than about 0.8 keV (see Fig. [3]). 

3.3. Tests of the extended X-ray structure 

Although the X-ray feature lies along the readout di- 
rection of the CCD, extensive testing supports our con- 
clusion that it is not an artifact but is instead repre- 
sentative of the actual X-ray morphology of RS Oph. 
Any excess charge from a readout artifact would ap- 
pear as a string of spurious X-ray events which, for 
bright sources (Chandra POG v. 11), are not removed 
by standard grade filtering and fall along rows of pix- 
els with constant y-coordinate in the CCD pixels array 
(CHIPY). Readout streaks are therefore always symmet- 
ric and might appear along pixels with the same CHIPY 
value as the source. Fig. [1] shows that the detected X- 
ray feature is not symmetric with respect to the central 
source. Furthermore, any readout artifact in the region 
of the CCD where we detect the feature would account 
for only a few counts, whereas we detect ~ 100 counts in 
the ~ \" .2 X-ray feature; we therefore conclude that it 
is not due to a readout artifact. 

Charge Transfer Inefficiency (CTI) in the S3 chip 
changes the overall instrument performance in various 
ways. First, since some charge is trapped, the charge 
read out is less than the amount of charge deposited. 
Second, CTI causes a migration of the grades assigned 
to each event due to photons that were captured by the 
traps and gradually re- emitted in the following frames. 
This second effect can change grades from that associ- 
ated with a single pixel events (flight grade = 0) to those 
associated with a vertical up- split (flight grades = 2, 66), 
which are oriented parallel to the readout direction (R. 
Edgars, private communication). Since no tested CTI- 
correction algorithm is available for observations taken 
in subarray mode, we performed several tests to indi- 
cate the degree to which our observation was affected by 
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Fig. 1. — Comparison of Chandra images of RS Oph with various levels of image processing, (a) Image obtained from pipeline products 
without further processing; (b) Image created after applying Subpixel Event Repositioning (SER, see text) to the data; (c) Resulting image 
after performing maximum likelihood deconv oluti on using a simulated PSF; (d) Image created after applying both SER and maximum 
likelihood PSF deconvolution as described in ^3.21 A comparison of panels (a) and (b) confirms that SER does not introduce any spurious 
asymmetry into the image. A comparison of panels (c) and (d) shows that the extended structure was clearly visible in the deconvolved 
image, even before the application of SER. The innermost contour encloses 4095, 3739, 4955, and 5140 counts in panels (a) through (d), 
respectively. Between the innermost and middle contours are 1252, 1589, 649, and 489 counts for panels (a) through (d), respectively. 
Between the middle and outermost contours are 299, 326, 66, and 64 counts for panels (a) through (d), respectively. North is up and east 
is to the left. 



CTI 6 . 

The S3 back-illuminated (BI) chip has substantial CTI 
in the serial array. If severe, this can create an extended 
feature in the direction along which the charges move 
out of the frame store region. In our observation, such 
an effect would correspond to a feature aligned in a direc- 
tion perpendicular to the one we detected. On the other 
hand, the detected feature could be caused by parallel 
CTI. To investigate the presence of parallel CTI, we ap- 
plied our image processing procedure to an observation 
from the Chandra archive that was taken with the same 
instrumental configuration (see § [2]) and has roughly the 
same count rate as our observation. The observation we 
used was of HD1 13703 B (ObsID 626, exposure time = 
12.7 ks). Trapped charges are unlikely to be re-emitted 
as grade-0 events. Since the narrow extended structure 
from RS Oph appeared strongly when we created an im- 
age using only grade events, we examined an image 
of HD1 13703 B created using only grade events. The 
grade-0 image of HD113703 B did not show any extended 
structure along the read out direction. This last test sug- 
gest that CTI effects were absent. 

We also carefully tested whether the image processing 
that we performed could have introduced the detected 

6 The CTI correction applied by the standard pipeline process 
was removed before performing any tests. 



feature as an artifact. SER is a event-by-event proce- 
dure that should not reposition the events in a prefer- 
ential direction. Nevertheless, to confirm that SER did 
not introduce an artifact, we applied the deconvolution 
procedure to data with and without SER applied (see 
Fig. [Tj). No significant differences were found between 
the SER and non-SER deconvolved images. We there- 
fore conclude that the narrow X-ray feature was not a 
by-product of applying SER. If the synthetic PSF was 
not symmetric, the deconvolution could have introduced 
artifacts in the resulting image. To test for this possi- 
bility, we applied the deconvolution with various rota- 
tions and flips of the PSF with no significant changes in 
the orientation of the extended feature. Finally, as we 
discuss below, features at roughly the same P.A. as the 
extended X-ray feature have been seen in radio and IR 
observations of RS Oph. 

4. DISCUSSION AND CONCLUSION 

Assuming, as suggested by our extensive testing, that 
the extended X-ray structure is not an artifact, the prop- 
erties of this feature indicate that it consists of a few 
tenths of a percent to a few percent of the ejecta mass, 
and that the outermost X-ray emitting material is mov- 
ing at a speed approximately half of the escape veloc- 
ity from the WD (the escape velocity is approximately 
12,600 km s" 1 for a 1.35 M WD). If the material com- 
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Fig. 2. — Profiles of the intensity along, and perpendicular to, the 
narrow X-ray feature. The solid curve shows the intensity profile 
along the narrow feature while the dotted line shows the profile 
in the direction perpendicular to the feature. The length of the 
extended feature, as measured from the centroid of the PSF to 
the location where the intensity falls below 10 counts per 0" .125 
spatial bin (marked with a dot-dashed vertical line), is 1.2±0.3 
arcsec. Error bars represent one standard deviation. The error in 
radial distance was estimated to be the best achievable SER pixel 
size, i.e., /x .3. 
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Fig. 3. — Images (to which both SER and maximum likelihood 
deconvolution have been applied) in 4 energy bands containing 
roughly the same number of counts (1000). Lower panel: Chandra 
spectrum of RS Oph (and overlaid best-fit model) with dashed 
vertical lines indicating the 4 energy ranges used to generate the 4 
images in the upper panels. The asymmetric X-ray feature is most 
pronounced in the soft band (E < 0.77 keV). 



prising the narrow X-ray feature was ejected at the time 
of the nova explosion, the length of the X-ray structure 
suggests an average expansion rate in the plane of the sky 
of ~2.3 mas d _1 . This rate is similar to the early expan- 
sion r ate of ~ 2 mas d _1 for the radio jet ( Sokol oski et al.1 
2008). The rate of expansion of the X-ray feature cor- 
responds to velocities in the plane of the sky of about 
- 6,300 ± 1600 km s' 1 (d/1.6 kpc). Since the X-ray 



emitting material is unlikely to be moving exactly in the 
plane of the sky, the true velocity could be higher. 

A strong shock moving at such a high speed would be 
expected to produce hard X-ray emission (with kT ~ 40 
keV in the strong-shock case); however we find only soft 
X-ray emission from the extended feature ( §3.2)) . The 
discrepancy between the expected temperature in the 
strong shock case and the observed photon energy is also 
prese nt in the diffuse X-ray e mission from planetary neb- 
ulae ( So ker fc Kas tner 200 31) and the two symbiotic stars 
with X-ray jets - CH Cyg (IGallowav fc Sokoloskill2004f ) 
and R Aqr (jKellogg et alJl2QQ7h , 

The narrow X-ray feature is most extended in the 0.3- 
0.8 keV energy band and contains approximately 40 pho- 
tons in this band. A crude estimation of its mass can 
be made if we assume that the emission is due to an 
absorbed, optically thin ther mal plasma with kT=0.5 
keV and n H =1.8xl0 21 cm" 2 (|Bode et al.ll2006f ). We use 
PIMMS to obtain the observed emission measure of such 
a plasma and, assuming that the narrow X-ray feature 
has a radius of - 0".9 and length of 1".2 (Fig. [2j), find 
that the density is ~40 cm -3 and the mass of the ex- 
tended X-ray emitting material is ~10 -8 M (assuming 
that there is one proton for every free electron). For 
comparison, estimates of the amount of material ejected 
immediately after the outburst range from 2jj_0_Z to 
-10- 6 M© (jSokoloski et al.1 12005 lHachisu et al.ll2007h . 

The orientation of the narrow X-ray feature is less 
consistent with that of the radio jets (lO'Brien et al.l 
120061 : iRupen et af]l2008l : iSokoloski etaIll2QQ8h than with 
other features detected at radio and infrared (NIR) wave- 
lengths. Early radio observations showed that the ring 
of synchrotron-emitting plasma associated with the ex- 
panding blast wave had a bright spot at a P.A. ^120° 
(jO'Brien et "all 120061 : IRupen et afll2QQ8f ). If one were to 
superpose the X-ray and radio images and draw a line 
along the extended X-ray structure, it would intersect 
the bright spots on the radio ring. Moreover, the radio 
ring was actually more of a C-shape, with the opening 
in the ring at the same P.A. as the X-ray feature. NIR 
interferometry uncovered extended emission with a P.A. 
of 120 o /300 ° from day 5.5 to a few months after the nova 
explosions ([C hesneau et all 1 2 0071 : [Lane et al.l [2007). and 
iChesneau et al.l (j2007l ) suggested a possible connection 
between the radio and IR morphologies. Our detection 
of an X-ray feature extending in the northwest direc- 
tion, like features in the radio and NIR, suggests that 
the structures along this direction could all be related 
(see Fig. [p. 

Although the X-ray structure has some similarities 
to features predicted by models of the explosion, many 
properties of the observed narrow X-ray feature are unex- 
plained, and the true origin of this extended X-ray emit- 
ting gas remains unclear. Computational modeling sug- 
gests that the blast-wave expansion should be strongly 
asymm etric. For example, in the 3-dimensional simula- 
tions of lWalder et al.l (|2008l ). the blastwave moves into a 
highly inhomogeneous pre-outburst density distribution 
and forms an elongated, bipolar structure perpendicular 
to the plane of the orbit. We do not see any extended 
X-ray emission in th is direction. In contrast, models by 
lOrlando et al.l (|2009l ) predict a one-sided X-ray structure 
- as we have found - due to the reflection of the blastwave 
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off of the red giant. Unfortunately, however, the position 
of the red giant at the time of t he nova explo sion (accord- 
ing to the ephemeris of iBrandi et aLlfeOOQl ) would have 
led to an X-ray extension in the north-south direction, 
which is inconsistent with our detection of extended X- 
ray emission at a P.A. of 300 degrees. Thus, our results 
suggest that either the orbital solution for this system 
needs to be re-examined or that the explosion somehow 
ejected X-ray emitting plasma in a direction that was 
neither perpendicular to the plane of the orbit nor along 
the line joining the WD and the red giant. 
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Fig. 4. — Compariso n of RS Oph im ages in radio, infrared and X-rays. Left: an image of radio synchrotron-emitting plasma with NIR 
contours overlaid (Chesneau et al. 2007). Right: the X-ray image obtained after SER and PSF-deconvolution processing (same as panel 
(d) in Fig. [T}. The X-ray feature extends along the same direction as the NIR and radio features, suggesting that they could all be related. 
Note that the two panels show very different scales. 



